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Abstract 

We study properties of the additive binary Markov chain with short and long- 
range correlations. A new approach is suggested that allows one to express global 
statistical properties of a binary chain in terms of the so-called memory function. 
The latter is directly connected with the pair correlator of a chain via the integral 
equation that is analyzed in great detail. To elucidate the relation between the 
memory function and pair correlator, some specific cases were considered that may 
have important applications in different fields. 

Key words: binary Markov chain, memory function, correlated disorder 
PACS: 05.40.2a, 02.50.Ga, 87.10.1e 



* F.M.Izrailev 

Email address: izrailevOsirio . ifuap .buap .mx (F. M. Izrailev). 



Preprint submitted to Elsevier Science 6 February 2008 



1 Introduction 



Much attention has been recently paid to the anomalous transport in ID 
systems with correlated disorder. In contrast with the well studied case of a 
white-noise potential, the problem of propagation of either electron or electro- 
magnetic waves through samples with random potentials with short or long- 
range correlations, is a big theoretical challenge. 

The fundamental significance of this problem is due to recent exciting findings 
revising a commonly accepted belief that any randomness in ID structures re- 
sults in the Anderson localization. In particular, it was shown [1,2,3,4] that 
specific long-range correlations in ID random potentials give rise to an emer- 
gence of frequency windows of a very high and very low transparency of waves 
through such potentials. It was also found that the positions and widths of 
these windows can be controlled by the form of the binary correlator of a 
scattering potential. A quite simple method was suggested for constructing 
random potentials that result in any predefined frequency window of an al- 
most perfect transparency. 

From the experimental viewpoint, many of the results may have a strong im- 
pact for the creation of a new class of electron nanodevices, optic fibers, acous- 
tic and electromagnetic waveguides with selective transport properties. The 
predictions of the theory have been verified experimentally [5,6] by studying 
transport properties of a single-mode electromagnetic waveguide with point- 
like scatterers. The agreement between experimental results and theoretical 
predictions was found to be unexpectedly good. In particular, the predicted 
windows of a complete reflection alternated by those of a good transparency, 
were clearly observed, in spite of strong experimental imperfections. 

A further development of the theory of correlated disorder is due to the appli- 
cation to the surface scattering of electromagnetic waves (or electrons). The 
problem of a wave propagation (both classical and quantum) through guiding 
systems with corrugated surfaces has long history and still remains a hot topic 
in the literature. For single-mode waveguides with the surface scattering the 
problem is equivalent to the ID bulk scattering. For this reason, the meth- 
ods and results obtained for the latter case can be directly applied for the 
waveguides [7,8,9]. Specifically it was shown, both analytically and by direct 
numerical simulations, that single-mode waveguides with a desired selective 
transport can be fabricated by a proper construction of long-range-correlated 
random surfaces. 

The key ingredient of the theory of correlated disorder is the binary correlator 
of a random potential. As was shown, for a weak potential this correlator fully 
determines the transmission/reflection of waves. The algorithm proposed in 
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Refs. [1,2,4] generates a statistical ensemble of random potentials having the 
same binary correlator. Note that in this method the values of the potentials 
are determined by a continuous distribution function. On the other hand, in 
some applications it is more appropriate to generate a correlated potential 
from a finite number of elements. An important example of such potentials is 
a sequence of nucleotides in a DNA molecule, where only four elements deter- 
mine a potential. In this paper we address the mathematical problem of the 
construction of a dichotomous sequence with a prescribed type of correlations. 
A priory, it is not clear whether the random sequence of only two elements 
may have a given binary correlator. To shed light on this problem we con- 
centrate our attention on statistical properties of a binary Markov chain and 
their relation to the binary correlator. Our main interest in this study is the 
so-called memory function that was recently analyzed in other applications 
(see Refs. [10-11] and references therein). We hope that our results may help 
to extend the theory of correlated disorder to binary potentials. 

The structure of this paper is as follows. In next Section 2 we introduce an 
additive binary Markov chain via the memory function. We show that the 
memory function can be related to the pair (binary) correlator of a chain. 
In Section 3 the integral equation for the binary correlator is developed and 
discussed in details. Complimentary, in Section 4 the equation for the memory 
function is derived in its general form. Next Section 5 is devoted to few other 
forms of the integral equations for the memory function, that can be more suit- 
able in specific cases. As the simplest case, in Section 6 we consider the Markov 
chain with a white-noise. Then, in Section 7 a bit more complicated situation is 
discussed, namely, the binary chain with short-range correlations of a general 
type. Another interesting case of an exponentially decreasing binary correlator 
is the subject of Section 8. Unexpectedly, in this case the memory function 
has a very simple, however, a quite specific form. Another modification of an 
exponential pair correlator, namely, modulated by the cosine-function, is ana- 
lyzed in Section 9. In next Section 10, we discuss a general case for which the 
Fourier transform of the pair correlator has single pole in the complex plane 
of its argument. Finally, in Section 11 we concentrate our attention on a very 
important and quite complicated case of the step-wise power spectrum. This 
case is of great importance in view of possible applications of our theory to 
the transport in physical or DNA systems described by the binary potentials. 



2 Additive Binary Markov Chain 

In what follows, we consider a homogeneous random binary sequence of sym- 
bols, 

e(n) = {e ,e 1 }, n e Z = -2, -1, 0, 1, 2, ... (1) 
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To specify an N-step Markov chain we introduce the conditional probability 
function, 

P(e(n) = £ ,i I e(n - 1), e(n - 2), e(n - N + 1), e(n - N)). (2) 



It is a probability of appearance of one of two symbols, e{n) = e or e{n) = ei, 
after a given sequence T N;7l ("word") of length N, 

T N , n = e(n - 1), e(n - 2), . . . , e(n - N + 1), e(n - N). (3) 



The additive Markov chain is defined by the conditional probability function 
of the form, 



P(e(n) = £ ,i I T N , n ) = p ,i + E F( r f in r) (4) 



r=l 



£0,1 — £l,0 



Here p and are the parameters that imply, respectively, the probability 
of occurring symbol e and e 1 in the whole sequence. Evidently, p + pi = 1. 
For a stationary random sequence these values po and p\ are, in fact, the 
only values of the probability density of symbols e(n). Note that the density 
does not depend on n, due to statistical homogeneity. Here and below the bar 
stands for the average along a sequence. In particular, 



M 



£{n) = a&oc 2MTT E e(n) = eoPd + eiPi. (5) 

1 n=-M 



One should stress that in our case this average is equivalent to the ensemble 
average. 

The function F(r) describes the effect of correlations between the n— th symbol 
e(n) and iV previous symbols e(n — r) with r = 1, . . . , N. In the following we 
refer to F(r) as to the memory function, since it may be associated with 
a memory about statistical properties of a sequence, necessary to generate 
random sequences with given characteristics. Due to the evident condition 
< -P(-l-) < 1, this function obeys the restriction, 

N 

J2\ F (r)\ < min(po,Pi)/max(p ,Pi). (6) 

r=l 



Our main interest is in the relevance of the memory function F(r) to such a 
statistical characteristic of random sequences as the binary correlation func- 
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tion, 

2 

C(r) = e{n)e{n + r) — e(n) . 



(7) 



The statistical homogeneity of the random sequence e(n) provides the cor- 
relator C(r) be dependent only on the distance r between two points n and 
n + r. Moreover, by definition (7) C(r) is an even function of this distance, 
C(—r) = C(r). From Eqs. (7) and (5), it follows that the variance C(0) reads 

C(0) = e 2 (n) - e(n) = (ei - e ) 2 p p!, with Po + Pl = I. (8) 

As was shown in Refs. [10,11], for additive Markov chains the relation between 
the memory function F(r) and correlation function C(r) has the form, 

N 

C(r) = F(r')C(r -r'), r > 1. (9) 

r'=l 

This equation provides a possibility to obtain the memory function F(r) from 
the prescribed correlator C(r) and to construct effectively the Markov chain 
with the conditional probability P{. \ .) according to Eq. (4). 

For simplicity, below we use the binary correlator K(r) normalized to the 
unity at r = 0, 

K(r) = C(r)/C(0). (10) 

Also, we assume, without the loss of generality, that the length of considered 
sequences is infinite. As a result, the starting equation (9) gets the form, 

oo 

K(r) = J2 F(r')K(r -r'), r > 1; (11) 

r'=l 

K(-r)=K(r), K(0) = 1. (12) 

Note that here the memory function F(r) is defined for r > 1, while the 
correlator K(r) being an even function, is determined for all values of r. To 
avoid mathematical problems related to this fact, we assume that the memory 
function F(r) can be analytically continued to negative values of r, namely, 
F(—r) = F(r). As one can see, we have also to specify the value of F(r) for 
r = 0. In what follows, for simplicity we assume that F(0) = 0. As a result, 
we arrive at main equations that shall be analyzed below, 
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K(r)= J2 K(r-r')F(r')-J2 K ( r + r ') F ( r ')i r > ^ ( 13 ) 

r'=— oo r'=l 

X(-r)=X(r), tf(0) = l; (14) 

F(-r)=F(r), F(0) = 0. (15) 



The above equations can be considered as the starting point for two prob- 
lems. The first one is to obtain the memory function F(r) by making use of 
Eq. (13) in which the correlator K(r) is known. The second problem is the 
complimentary one, namely, to find the correlator K(r) for a given memory 
function F(r). 



3 Equation for Binary Correlator 



It is convenient to write the Fourier representation of the main equations (13). 
For this, we introduce the Fourier transform K,(k) of the binary correlator 
K(r), known as the randomness power spectrum, 



1 *" 

K{r) = - ( dk)C(k)cos(kr), (16) 

7T J 



IC(k) = 1 + 2j2K(r)cos(kr); (17) 



r=l 



K{-k)=K{k), }C(k + 2n) = K{k). (18) 



Since the correlator K(r) is a real and even function of the coordinate r, its 
Fourier transform (17) is an even, real and non-negative function of the wave 
number k. The condition Kir = 0) = 1 results in the following normalization 
relation for the power spectrum tC(k), 

— I dklC(k) = - f dkK{k) = 1 . (19) 

2,71 J TV J 



Analogously, the expressions for the Fourier transforms of the memory function 
F(r) have the form, 
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1 n 

F(r) = - [ dkF(k)cos(kr), (20) 

71 J 



F(k) = 2^F(r)cos(Ar); (21) 



r=l 



F(-k) = F(k), T{k + 2n)=T{k). (22) 



Similarly, the Fourier transform !F{k) of the memory function F(r) is also 
even and real function of the wave number k. Due to the imposed condition 
F(r = 0) = 0, the normalization condition for J-'(k) reads as, 

- / dkT{k) = - [ dkF(k) = 0. (23) 

'.71 J 71 J 



2tt 

— 7T 



Now, we shall derive the Fourier transform of relation (13). First, we represent 
this equation in the following from, 



K{k) = 1 + Yl exp(ikr)K(r - r')F(r') - ]T K(r)F(r) 

r,r'=-oo r=— oo 

oo 

-2 cos(kr)K(r + r')F(r'). (24) 

r,r'=l 

Now we substitute the memory function F(r') due to Eq. (20), together with 
the use of Eqs. (17). In this way we arrive at the relation, 



K,{k)F(k) - 2 jr cos(kr)K(r + r')F(r') = K,{k) + A - 1, (25) 

r,r'=l 

where the constant A is defined by 

1 r 

A= ^ K(r-)F(r-) = — J dkK{k)F{k). (26) 



Into the second term of the left-hand side of Eq. (25) we substitute the cor- 
relator K{r + r') in the form of the Fourier integral (16). As a result, we get 
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the integral equation for the power spectrum lC(k) of the binary correlator, 

7T 

[T(k) - l]K{k) - — I dk' T(k, k')K,(k') = A - 1. (27) 
2ir J 

Here the integral kernel T(k, k') is described as follows, 



oo 



T(k, k') = 2j2 cos(kr) exp(-ik'r) x J2 F ( r ') exp(-ifcV) (28) 

r=l r'=l 

1 l r 

= o / dk"F{k")Q{k", k'), k' = lim (k' - ie), 

2 7T J e— >+0 



and the Q-functions are defined below by Eq. (34). 

Due to a complex nature of Q-functions, the kernel T(k, k') is also a complex 
function. Then, in accordance with the definition (28), we have, 



7T 7T 

— [ dkT(k,k f ) = - fdkT(k,k') = 0, (29) 

27T J 7T J 

-7T 

7T TT 

- / dfc' T(Jfe, k') = - ( dk' T(k, k') = 0. (30) 



2vr 

— 7T 



When integrating Eq. (27) over the wave number k within the interval (— tt, it) 
or (0,7r), its integral term vanishes due to Eq.(29), while the equation itself 
turns into identity A — 1 = A — 1. This fact implies that the relation (26) 
automatically satisfies the solution of Eq. (27), independently on the value 
of A. Therefore, the constant A is determined solely by the normalization 
condition K(0) = 1 or, the same, by Eq. (19) for the power spectrum K{k). 
Another form of Eq.(27) will be given below by Eq. (53). 



4 Equation for Memory Function 



Let us now substitute the memory function (20) into the left-hand side of 
Eq. (25). In this way we obtain the integral equation for the Fourier transform 
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J-{k) of the memory function, 

K{k)F{h) - - I dk'E(k, k')F(k') = K,{k) + A - 1. 

TT J 



(31) 



In this equation the kernel E(k, k') of the integral operator is described by the 
expressions 



E(k,k')=2 J2 cos(kr)K(r + r') cos(k'r') 



r,r'=l 



1 

- J dk'%(k")Q(k,k")Q(k',k"), 



where the Q-functions are defined as 



(32) 
(33) 



Q(k, k") = -2j2 cos(kr) exp(-ik"i 



= l + i 



r=l 

sin k" 



k"= \iui(k" -it). 
cos k — cos k" 



(34) 



(35) 



In spite of a complex nature of the Q-functions, the kernel E.(k,k') is real, 
even and symmetrical function of both arguments k and k', 



E{-k, k') = E(k, —k') = E(k, k') = E(k', k) 



(36) 



One can show that the mean value of the function Q(k, k") vanishes, 



i i 

— / dkQ(k,k") = - ( dkQ(k,k") = 0. 

-TV 



(37) 



Indeed, integrating Eq.(35) over k, one can obtain, 



sin k"dk 



ITT 



2k c J cos k — cos k" k c 

—k c 



Q(k c -\k"\) + 2 f;^&^sin(A"/ 



fe(k c -\k"\) + Un 



sin[(A; c + k")/2] 



sin[(A; c - k")/2] 



i k c r 

0<k c <7r; \k"\ < 7i. (38) 



9 



The first expression for the integral is obtained by a direct integration of the 
sum in Eq.(34). In the derivation of the second expression we have used the 
fact that the Heaviside unit-step function Q(x), for which Q(x < 0) = and 
Q(x > 0) = 1, arises due to the contribution of simple poles k = ±(k" — ie), 
e — > +0. Clearly, this contribution emerges only when the poles are within the 
integration interval, i.e. if \k"\ < k c . The principal value (P. V.) of the integral 
in Eq.(38) has the following form, 



p.v. f sinadx =- P.V. [ 

J cos x — cos a 2 J 
o 

sin[(fc c + a)/2] 



dx 



cot 



x + a 



cot 



x — a 



In 



sin[(/c c — a)/2] 



(39) 



One can see that the integral (39) vanishes when k c = ir, and the Q-function 
property (37) is evident. 

Note that Eq.(37) gives rise to the normalization relations for the kernel 
E(k,k'), 



7T 7T 

— [ dk Z(k, k') — — f dk E(k, k') — 0, (40) 

2,71 J 71 J 

-7T 
1 ^ 1 ^ 

— / dk' E(k, k') — — f dk'E(k, k') = 0. (41) 

Z7l J 71 J 

-7T 

As a result, the integration of Eq. (31) over the wave number k within the 
interval (— n, n) or (0, tt) leads to the identity A = A. This fact implies that 
the relation (26) automatically satisfies the solution of Eq. (31). Therefore, 
the constant A is determined by the initial condition for the memory function 
F(r), 

7T TV 

F(0) = — [ dkT{k) = - [ dkT{k) = 0. (42) 

2,71 J 71 J 



5 Other Forms of Equation for Memory Function 



In some cases, instead of Eq. (31) it is more convenient to use other forms of 
the integral equation for the Fourier transform T{k) of the memory function. 
Specifically, the expression (32) for the integral kernel E(k, k') turns out to be 
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ineffective for some kinds of the correlator K(r). In this case one should use 
Eq. (33) and, therefore, deal with complex Q-functions (34), taking explicitly 
into account their pole contributions. In what follows, we assume that K,(k) 
and J-{k) are even, periodic and sectionally continuous functions on real axis 
k, see Eqs. (18) and (22). 

In accordance with Eq. (35) the product of two Q-functions can be written as 



Q(k, k")Q{k\ k") = 1 + i f nfe " - (43) 

cos k — cos k" 

sin k" sin 2 k" 



cos k' — cos k" (cos k — cos k") (cos k' — cos k") 

After the substitution of Eq. (43) into the definition (33) for the kernel E(k, k'), 
one can obtain, 

2 E(k, k') = l- K(k) - K(k') - 2E {k, k'). (44) 



Here the first term arises due to the unit in Eq. (43) and normalization condi- 
tion (19). The second and third terms are direct consequences of the following 
relation, 



— ( dk"K,(k") ^ = %K(k), \k\ < n. (45) 

2tt J cos A; - cos k" w 1 1 ~ v 1 

This relation can be obtained from a more general one, 



( dk"K{k") * mk " - = ^Q(k c - \k\))C(k), 

J y ' cos A; - cos k" K y 1 u y h 



2 



0<k c <n, \k\ < TT. (46) 



The integral in Eq.(46) is determined by simple poles k" — ie = ±|/c|, e — > 
+0, if the poles are within the integration interval, i.e. when \k\ < k c . The 
corresponding principal value of the integral vanishes because the integrand 
is an odd function. Thus, when k c = tt, Eq. (46) reduces to Eq. (45). 

One can see that in Eq. (44) a new integral kernel 

E (k, k') = — J dk"JC(k") (cog - - ^ (cog - - W) (47) 
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is introduced. Similar to the function E(k,k'), the new kernel Eo(k,k') is a 
real, even and symmetrical function with respect to both arguments k and k', 

E (-k, k') = E (k, —k') = E (k, k') = E (k', k). (48) 



In line with Eq. (38) taken at k c = n and Eq. (45), this new kernel E (k, k 1 ) 
satisfies to the following integral properties, 



\ *W , (49) 

-±K(k). (50) 

These relations, together with Eq. (44) and normalization condition (19), pro- 
vide the properties (40), (41) of the old kernel E(k, k'). 

Finally, with the use of Eq. (44), the normalization condition (42) for the 
memory function, and the definition (26) for the constant A, we arrive at the 
new equivalent integral equation for the memory function in /c-representation, 

7T 

lC(k)F(k) + - I dk' E (k, k')F(k') = K{k) + -A-1. (51) 

TT J 2 




dkE (k, k') = - 
dk'E (k,k') = - 



As one should expect, when integrating Eq. (51) over the wave number k 
within the interval (— 7r, 7r) or (0, tt), its integral term gives the value — A/2 due 
to the equality (49), thus, reducing the equation into the identity, A/2 = A/2. 

Now let us derive one more form of the integral equation. To this end, one 
should substitute the expression (47) into the integral term of Eq. (51), with 
a further change of the order of integration over k and k", 



- / dk'Z (k,k')F(k') 
ix J 

IT 

J dk"K(k" 



1 

47T 



sin k" 1 

x 



cos k — cos k" tt 



dk' F(k') 



sin k" 



cos k' — cos k" 



sin K 



k" 



cos k — cos k" 



X 



1 

iF(k") + - P.V. ( dk'F(k'] 

7T J 



sin K 



k" 



cos k' — cos k" - 
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1 1 J 

= --K{k)T{k) + — j dk'%{k") 



sin k" 



cos k — cos k" 



x- P.V. f dk' Hk') y^" ,„ ■ (52) 
7T J cos k' — cos k" 





Here we have employed the ideas used to evaluate the integral (38), and the 
formula (45). After we substitute the result (52) into Eq. (51), the integral 
equation reads 



7T 

mm + 7^ dk'%(k") smfc - (53) 

2tt 2 J cos k — cos k" 

— 7T 

x P.V. ( dk'F(k') ^ = 2K{k) + A-2. 

J cos k — cos k" 



As above, the symbol " P. V." denotes a principal value of the corresponding 
integral. Since we cannot change the integration order in Eq. (53), its applica- 
tion is suitable if the power spectrum JC(k) is considered as unknown while the 
memory function F{k) is predefined. When integrating Eq.(53) over k within 
the interval (— ir, it) or (0, it), its integral term vanishes due to Eq. (38), while 
the equation itself turns into the identity A = A. 

One should, however, emphasize that all the integral equations (31), (51), and 
(53) cannot be solved in general case of any form of the additive-Markov-chain 
power spectrum K{k). This is because of a quite complicated structure of the 
integral kernels. For this reason, below we shall analyze some particular cases 
when the solution can be obtained in a relatively simple way. 



6 White-Noise Disorder 



The white-noise Markov chain is specified by the binary correlator and corre- 
sponding power spectrum as follows, 



K wn {r) = 5 rfi , JC wn (k) = 1. (54) 

In order to obtain the memory function F wn (r) of the additive chain, one 
can use the equation (11), or (13). Due to the restriction r > 1, one should 
substitute zero in their left-hand sides, that gives rise to the expected result, 

F wn {r) = 0. (55) 
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Now, let us see how the same result follows from the integral equations (31), 
(51) and (53). One should take into account that for a white noise with 
fcwn(k) — 1 the definition (26) for A and the normalization condition (42) 
yield A = F(0) = 0. Therefore, the right-hand side of the equations vanishes 
and they become homogeneous, with a trivial solution, 

r wn {k) = 0. (56) 



In such a way, we again come to the result (55). 

Let us write down the integral kernels of the equations with the white-noise 
power spectrum K. wn (k) = 1. First, for Eq. (31) we have, 



Z wn (k,k') = 2 cos(kr)5 r+r/fi cos(k'r') 

r,r'=l 

= — J dk"Q(k, k")Q(k', k") = . (57) 

In this expression the sum is evidently equal to zero. To make sure explicitly 
that the integral also vanishes, one should employ Eqs. (43), (45) and the 
relation 



27r J (cos k — cos k") (cos k' — cos k") 



As for Eq. (51), taking into account Eq. (58), we get 

SHM'H-1/2- (59) 



7 Short-Range Correlations 



The simplest non-trivial case for which one can obtain an asymptotically exact 
solution of Eq. (31), is the additive Markov chain with short range correlations. 
In this case the correlator K(r) is assumed to be a rapidly decreasing function 
of the index r, with a "very short" scale R c > of its decrease. It is evident that 
for short-range correlations the power spectrum K,(k) is a very smooth function 
of the wave number k with a "long" scale k c = R~ l > of decrease. One should 
take into account a discrete nature of the argument r of the correlator K(r), 
or, the same, the periodicity (18) of its power spectrum JC(k). Then, it becomes 
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clear that depending on a specific form of K(r) (or K,(k)), the definition of 
short range correlations assumes one of two conditions, 



R c = k c 1 < 1, or 7r - k c < 7r. (60) 



In this case the kernel S(fc, fc') of the integral operator, see Eq. (32), is mainly 
contributed by values r = r' = 1. Consequently, one can get the following 
estimate for S(/c, fc'), 

S(/c, k') « 2tf(2) cos fc cos fc' < X(0) = 1. (61) 



Therefore, the integral term in Eq. (31) can be neglected. As a result, the 
integral equation transforms into an algebraic one with the solution, 

F{k) = 1 - (1 - A)Kr\k). (62) 



By making use of Eqs. (20) and (62), the memory function in the coordinate 
representation reads 

F(r) = 5 r , -(l-A)I(r), (63) 



with the integral I(r) defined by 

— 7T 



Applying the initial condition (42), one can easily get 

1-A = r\0). (65) 



As a result, in the case of short-range correlations we arrive at a very simple 
expression for the memory function F(r), 

F(r) = 5 r , - /(r)//(0). (66) 



Evidently, this expression should have a crossover to Eq. (55) for a white-noise 
disorder. Indeed, if K.{k) = 1, one easily gets I(r) = 5 r> o and F(r) = 0. 
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8 Exponential Correlations 



Clearly, the equation (31) can be easily solved if the kernel E(k,k'), see 
Eq. (32), is factorized by a product of two terms that separately depend on k 
and k! only. 

A particular important example, giving rise to such a form of the kernel 
S(/c, k'), is the Markov chain with an exponential correlator, 



*.<'> = «P(-*Jr|), W- J^Xe* - < 67 > 

Here k c = R^ 1 > determines the inverse correlation length. The explicit 
form (67) for the randomness power spectrum K,(k) can be found with the use 
of the relation, 

2 £oo8(Jbr) exp(-A; c r) = - 1, (68) 



that can be obtained by a direct summation of the geometric progression, after 
rewriting cos(kr) in the Euler form. 

In accordance with Eqs. (32), (67), and (68), the kernel of the integral operator 
S(fc, k') can be written in the factorized form, 

~ e (k,k') = ^{JC e (k)-l][JC e (k')-l]. (69) 



This form provides a possibility to find easily the solution of Eq. (31), 

Hk) = lf^O) [1 " 1 /A0)/C e (A;)] , (70) 

F{r) = TT^o) [6r >° ~ ' ( 71 ) 



with I(r) given by the definition (64). 

Taking into account the expression (67) for K,(k) when evaluating the integral 
in Eq.(64), one can get 

I(r) = cothk c 5 r fi — (2 sinh k c )~ l 5\ r \ t i . (72) 
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Finally, we come to the expression for the memory function, 

F(r) = exp(-A; c )<J| r | ) i . (73) 



It is important that the additive binary Markov chain with an exponential 
correlator is nothing but a sequence with the one-step memory function. As 
one can see, the exponential nature of the correlator K(r) results in a one- 
step-form of the memory function F(r), while its variation scale determines 
the amplitude only. 

For the Markov chain of symbols {eq — 0, E\ — 1} with the relative number 
Pi of "1" in the chain, the conditional probability of the symbol "1" occurring 
after the symbol "1" is equal to 

P(e(i) = 1 | 1) = pi + exp(-& c )(l - Pl ). (74) 



Here we took into account that (s(i)) = Po x + pi x 1 = p 1 . In the limit of 
an anomalously small correlation length k~ x — > 0, the conditional probability 
function tends to constant value, P(e(n) — 1 | 1) — > p\ that is independent 
of previous symbols. This limit evidently corresponds to the uncorrelated se- 
quence. In the opposite case of long-range correlations, k^ 1 — > oo, the condi- 
tional probability function is almost equal to unity, P(e(n) — 1 | 1) — > 1. In 
this case we have a sequence with a maximal persistent diffusion. 

There is another important example for which the kernel of the operator 
H(/c, k') has a separable form. Namely, this is the Markov chain with an expo- 
nential alternating- sign correlator, 



sin h /r 

K(r) = (-l) r exp(-£; c |r|), £(*) = —- c —-. (75) 

w v ' v 1 u w cosh k c + cos k v 1 

For this case, we obtain the following integral I{r), 

I {r) = coth k c 5 r fi + (2 sinh /c c ) _1 5| r M , (76) 



and, correspondingly, the negative memory function, 



F{r) = TTW) [5r '° " J ^/ J ^] = " e M-k c )6\r\A ■ (77) 

As is known [10], the negative memory function describes the anti-persistent 
correlated diffusion in the additive binary Markov chain. 
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To conclude, we would like to note that in the limit of short-range correlations 
with k c ^> 1, one gets 1(0) ~ 1. This leads the exact solutions (73) and (77) 
that coincide with the asymptotic expression (66). 



9 Exponential-Cosine Correlations 



In previous Section 8 we have revealed that the exponential binary corre- 
lator (67) results in the one-step memory function (73). In order to have a 
many-steps memory function, we consider the oscillating correlator with an 
exponential decrease of its amplitude, 



K ec (r) = exp(-A;o|r|) cos(cur) = X - [K e (r) + K* e (r)}; (78) 



^, (cosh k — coscjcosA;) sinhA; 

(cosh k cos uj — cos k) 2 + sinh 2 k Q sin 2 uj 

= l -[K e {k)+K* e {k)\. (80) 

Here the asterisk "*" stands for the complex conjugation. The last expres- 
sions (78) and (80) are given to relate the exp-cosine correlator and its power 
spectrum to those (67) for purely exponential correlations. Based on this re- 
lationship, one should perform the following complex substitution, 

k c = ko + iuj. (81) 



In line with the definition (32) and representation (78), the kernel S(/c, k') of 
the integral operator is described by the sum 

Z ec (k,k') = ±[Z e (k,k') + Z* e (k,k% (82) 



Therefore, according to the factorized form (69) of the operator Z e (k,k'), 
the integral equation (31) gives rise to the following solution for the memory 
function, 



2 + A 2- A 1 B K, e (k) - K.* e (k) 
^ ec(A;) " ~2 2~ Kjk) + 2 2K ec (k) ' (83) 
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Fe C {r) = ^ 5 r , - ^ I(r) + | $(r) . (84) 

Here the real constant A is defined by the standard expression (26), while the 
imaginary constant B is defined by 

B = l/' a M ,. lt) , (85 ) 





For real function J(r) the standard definition (64) remains valid, and new 
imaginary function $(r) has the form, 



— 7T 

K e {k) — K* e {k) . (cos c<j — cosh /c cos k) sin a; 

2 (cosh A; cosc<j — cos A;) 2 + sinh 2 /c sin 2 cj 



In order to find the values of constants A and B, one should write down two 
equations. The first one is nothing but the normalization condition (42). The 
second equation can be obtained by substitution of the solution (83) into the 
definition (85). As a result, we come to the system 

[/(0) + 1]A + $(0)5 = 2[/(0) - 1] , $(0)A - CB = 2$(0) , (87) 
where the real constant C is given by the expressions, 



r _ 9 i 7 [JC e (fc)-jc;(fc)] 2 _ i j jc e (k)jc e (k) 

2^J dk TKJk) " 1 + 27 J dk K ec (k) 

— TV 

The solution of the system (87) is 



I(0)-l + $ 2 (0)/C 4$(0)/C 

J(0) + l + $ 2 (0)/C" " /(0) + l + $ 2 (0)/C 1 J 



Thus, the memory function F(r) of the additive Markov chain with exp-cosine 
correlations has the form, 



2 f $ 2 (0) $(0)$(r)i 

F »< r) = /(0 ) + i + gffl i |/(0) + T 1 '* - /(r) - J 4^ ii } ' (90) 
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One can see that for uj = the power spectrum is JC e (k) = fC*(k) = JC ec (k), 
so that we have $(r) = 0, and Eq. (90) coincides with Eq. (71). 



10 Power Spectrum with Simple Complex Poles 

It is evident that the Fourier integrals (16), (20) in the definitions of the 
binary correlator K(r) and memory function F(r) are completely determined 
by singularities of corresponding Fourier transforms JC(k) and J-{k) in the 
complex plane of k. Taking into account the periodicity conditions (18) and 
(22), it is clear that the considered Fourier integrals are determined by the 
singularities located within the complex strip (— tt, 7r). 

Let us consider the integral kernel E(k, k') represented by Eq.(33). It is clear 
that Q-functions in the integrand are analytical functions of the integration 
variable k" in the lower half-plane of k" . Indeed, they have periodically re- 
peated simple poles k" = ±|fc| + ie and k" — ±|fc'| + ie (e — > +0) solely in the 
upper half-plane and converge to zero when k" — > — ioo, due to the expression, 

Q(k, k")Q(k', k") « 4 cos k cos k' exp(-2y), 

for k" = —iy, y — > oo. (91) 

Therefore, the integral in Eq. (33) is essentially specified by the analytical 
properties of the power spectrum fC(k) in the lower half-plane. Let us evaluate 
the integral along the closed contour ABC DA with A = —n—i0, B = —n—ioo, 
C = 7t— too and D = n—i0 (see Fig. 1). The sum of integrals along the straight 
lines AB and CD is equal to zero since the integrand is a periodic function with 
the period 27r. Then, the integral along the infinitely far line BC also vanishes 
because the product of two last factors of the integrand vanishes exponentially 
fast at any infinitely far point in the lower half-plane, see Eq. (91). Thus, the 
integral (33) taken with the opposite sign, is completely determined by the 
singularities of JC(k) inside the considered contour ABCDA, or the same, 
within the lower part of the complex strip (— vr, n). 

Let the power spectrum JC(k) be analytical within the considered strip. In this 
case JC(k) is a constant and due to the normalization condition (19), it has 
to be one, JC(k) = 1. From our analysis it follows, S(fc, k') = 0, that is in an 
accordance with the previously obtained results for a white-noise disorder, see 
section 6. 

Now we are in a position to analyze the case when there is only one singularity 
within the strip. Specifically, we consider a simple pole at k" = —ik c , assuming 
that the function K,(k) can by represented in the vicinity of this pole by the 
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D Rek" 




Fig. 1. The contour of integration ABCD in the case of one singularity in the strip, 
form, 



K{k") 



k" + ik c 



(92) 



Here we take into account that the normalization condition (19) requires the 
residue be equal i. In this case one can get, 



E e (k, k') = ~Q(k, -ik c )Q(k', -ik c ) . 



(93) 



One can see that the above expression coincides with Eq. (69), thus reducing 
our problem to that considered in Section VIII. 



11 Step- Wise Power Spectrum 



11.1 General expressions 



As was mentioned in Section I, in the theory of correlated disorder [1,2,3,4], 
one of important problems is an additive binary Markov chain with a step-wise 
power (SWP) spectrum. The simplest case is given by the following expres- 
sions, 



l\(r) = i h)^P^. (94) 



k c r 
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K{k) = h+(l-h)^- [Q(k + k c ) - Q(k - k c )\ = h + (1 - h)A kc (k) , (95) 

0</i<l, < A; c < 7T, |A;|<7r. 

The restriction on the parameter h is a consequence of a non-negative nature 
of the power spectrum K,{k). Strictly speaking, the exact requirement is < 
h < n I '{it — k c ) . In our analysis we assume h be smaller than one, in order to 
allow k c be very small. 

In Eqs. (95) and (96) the unit-pulse function A kc {k) is defined by 



A kc (k) = ^Q(k c -\k\), 0<k c <n, \k\<7r, (97) 

IT IT 

— J dkA kc (k) = -j dkA kc (k) = 1, (98) 

-7T 

lim A kc (k) = lim ^-&(k c - \k\) = 27v5(k). (99) 

This definition and hence the power spectrum K,(k), together with its in- 
verse quantity /C _1 (&;), contains Q(x) as the Heaviside unit-step function, 
e(x < 0) = and Q(x > 0) = 1. The factor n/k c in Eqs. (95) and (97) 
provides the normalization condition (19). As the majority of unit-pulse func- 
tions, our A kc (k) is also a prelimit Dirac delta-function. This fact is displayed 
by Eq. (99). The Fourier transforms for A kc {k) are 



a /,\ sin(A; c r) ^^sin(k c r) . ,, nn \ 

&kM= E V ; exp(tfcr) = 1 + 2^2 I cos(kr); (100) 

r=-oo r=l ^ cT 

= J_ f dkA kc (k) exp(-iAx) = - f dkA kc (k) cos(kr). (101) 

k c T Z7T J 7T J 

-IT 

Note that Eq. (100) for A kc (k) coincides in form with the Fourier transforms 
(16) for K.{k) and K(r). This is not surprising since 



K(r) = sm(k c r)/(k c r), K{k) = A kc (k) for h = 0. (102) 

Therefore, one can expect that all formulas derived in this section have to be 
reduced to the corresponding general ones in the limit when h — 0, with the 
change, A kc (k) -> K(k). 
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One can see that in the case when k c — > 0, we have the sequence with long- 
range correlations, whereas if k c — > n the system is a short-range correlated 
and at k c = it it reduces to the white-noise (54). More generally, when k c 
diverges, (k c — > oo), the second term in Eq. (94) decreases with oscillations. At 
the same time, K{r) gets the white-noise form (54) every time when k c = nn 
(n = 1, 2, 3, . . .). We emphasize that Eqs. (95) - (96) are applicable only when 
k c < tt. Note that the crossover to the white-noise (54) can be also performed 
with h — > 1. 

In order to obtain the memory function F(r) for the Markov chain with the 
step-wise power spectrum (94), one should firstly derive the kernel S(/c, k') of 
the integral equation (31). For this, we substitute Kir + r') or JC(k") in the 
form (94) or (95) into the definitions (32) or (33), respectively. The first term 
with h vanishes due to the condition (57), therefore, we have 



-n <' 1 ~ h A -ST (i \ sin M r + r ')] / 7 / /\ 
^{k,k) = — - — 4 2^ cos(/cr) — j— — - cos(/c r ) 



1-h 1 



2 2tt 



1-h 1 



2 2n^^ 



ry=1 K(r + r') 

IT 

J dk"A kc (k")Q(k,k")Q(k',k") 
dk"Q{k, k")Q(k',k"). 



(103) 
(104) 

(105) 



Now we employ Eq. (43) for the product of two Q-functions to write explicitly 
some integrals. This yields, 



Z(k,k') 



1 - h 



1-A kc (k)-A kc (k') + I sw (k,k') 



(106) 



Here the new integral kernel reads, 



1 

I sw (k,k') = -— J dk"A kc {k") 



x 



sin 2 k" 



[cos k — cos k") (cos k' — cos k") 



~, (107) 



1 

~2k r 
1 - 



sin 2 k"dk" 



-k, 

1 1 



cos k — cos k") (cos k' — cos k") 

sinp c + k)/2\ 



(108) 



k r cos k — cos k' 



sin k In 



sin[(£; c - k)/2] 



— sin k' In 



sinp c + fc')/2] 



sin[(A; c - fc')/2] 



(109) 
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When k c = it, we have I sw (k,k') = 1, and the kernel E(k,k') vanishes. This 
fact is in an agreement with the condition (57). Actually, the same expression 
(106) for the kernel S(fc, k') can be achieved from Eq. (44) by taking explicitly 
into account the definition (47) and Eq. (58). 

The kernel I sw (k, k') is a real, even and symmetrical function of both variables 
k and k', 

Isw{ k-> k ) = I sw (k, k ) = I sw (k, k ) = I sw {k , k). (HO) 

In line with Eq. (38) taken at k c = tt and Eq. (46), this new kernel I sw (k, k') 
satisfies the following integral properties, 

- fdkI sw (k,k') = A kc (k'), (111) 

TT J 


- / dk'I sw (k,k') = A kc (k). (112) 

TT J 


These relations together with Eq. (106) and normalization condition (98) pro- 
vide the integral properties (40), (41) of the old kernel S(fc, k'). 

Now we substitute Eqs. (95) and (106) into Eq.(31). One should take into 
account that two first summands from E(k,k') do not contribute into the 
integral term of Eq. (31), due to the normalization condition (42). As a result, 
one can get the integral equation for the memory function, 



7T 

T — -jF(fc) + A kc (k)f(k) - — f dk'I sw (k, k')F(k>) 



= A kc (k) - 1+ l -B sw . (113) 

Here, instead of the old constant A, for the convenience we have introduced 
the new one, B SW) which can be associated with A according to Eq. (26), 



A =(l-h)B sw , (114) 
B sw = — j dkA kc (k)T(k) = — J dkT{k). (115) 



Note that Eq.(113) can be derived from Eq. (51). 
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Finally, let us integrate Eq. (113) over the wave number k within the interval 
(— 7T, 7r) or (0, 7r). Then, on its left-hand side the first term vanishes due to 
Eq. (42), second term turns into B sw in accordance with its definition (115), 
while the last integral term gives the value B sw /2 because of the property 
defined by Eqs.(lll) and (115). On the right-hand side, the first term com- 
pensates the unity. In such a manner, the equation is reduced to the identity 
B sw /2 = B sw /1. 

It should be, however, stressed that having Eq.(113) in its general form, it is 
still not clear how to derive the solution. 



11.2 SWP Spectrum: Short-Range Correlations 



A particular case is the SWP spectrum with short range correlations, which 
corresponds to the values of k c close to ir. In this case Eqs.(94) and (95) result 
in a particular case of Eq. (66) with the integral I(r) that can be calculated 
in accordance with its definition (64). Taking into account Eq. (96) one can 
derive, 



J(r) 



1 

h 



dr,0 — 



k c (l — h) sin(A; c r) 



7r — (tc — k c )h k c r 



(116) 



Therefore, the memory function has the form, 



F(r) 



k c (l - h) 



(tt - fc c )(l -h) + k c h 
for 7r — k c <C ir. 



sin(/c c r) 
k c r 



-5. 



r,0 



;ii7) 



As a result, within the first approximation in h — > we get, 



F(r) 



k r 



TT - k r 



sin (k c r) 
k r r 



(118) 



which is valid for k c h <C re — k c <^ ir. 



11.3 SWP Spectrum: Long-Range Correlations (Zero Approximation) 

Now let us consider the more interesting case of the SWP spectrum (95) with 
long-range correlations. In this case the additive binary Markov chain has 
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small wave number k c <C 1. Thus, within zero approximation in k c the exact 
expressions (95) for the binary correlator K(r) and its power spectrum K,{k) 
can be changed by the following asymptotics, 



K{r) = h5 rfi + {l-h), (119) 
K(k) = h + (l-h)2n5(k). (120) 

with \k\ < ir. Note that these expressions are in an agreement with Eqs. (97)- 
(99) and (100)-(101). 

Now we substitute K(r + r') or K{k") in the form of Eq. (119) or Eq. (120) 
into the definitions, respectively, (32) or (33) for the kernel E(k,k'). Taking 
into account that the first term with h vanishes due to the condition (57), we 
get the result 



1-/' - ■ • < - 121 > 



r=l " r'=l 

l-h 



E(k, k') = — — [2 cos(kr)] [2 ]T cos(A;V 

r'=l 

Q(k,0)Q(k',0) (122) 



= [1 - 2n5(k)} [1 - 2^(^)1 (123) 

= - [l - 2n5(k) - 2n5(k f ) + An 2 5(k)5(k f )] , (124) 

that should be compared with the exact expressions (103)-(105) and (106). In 
order to pass from Eqs. (121)- (122) to Eq. (123) and then to Eq. (124), one 
should recognize that from the definition (34) for the Q-function, it follows, 

00 

Q(k, 0) = -2j2 cos(At) = 1 - 2n5(k). (125) 

r=l 



Note that the expressions (121)-(124) and (125) fulfil the general relations 
(40)-(41). 

The integral equation (31) with Eq.(124) taken as the kernel S(fc, &'), gives 
rise to the equation, 



F{k) + 27c5(k)F(k) - n5(k)T(0) = 2n5(k) + - ^(0). (126) 



l-h v ' v ' v ' v ' v ' v/ 2 

Here we take into account that 

B sw = J^(0) for k c -> 0. (127) 
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All terms in the equation (126) correspond to those in the general equation 
(113). Thus, the solution for the memory function yields 



F(r) 



1 - h 
h 

1-h 

h 



1-^(0)][27T*(*)-1], 



1 - - F(0)\ [1 - 5 r ,o] , h^O. 



(128) 
(129) 



In principle, these relations would give a complete solution of the problem. 
However, a problem remains how to find the constant B sw = ^(O). Indeed, 
the solution (128)-(129) automatically satisfies to the normalization condition 
(42). On the other hand, the relation (128) violates if one takes k — 0. 



11. 4 SWP Spectrum: Long-Range Correlations (Another Approach) 



Our failure to derive the correct solution in previous subsection is evidently 
related to the fact that at the beginning we have substituted the rigorous 
expressions (94)-(95) by Eqs. (119)-(120). Strictly speaking, the latter are 
always invalid since the smaller k c the larger indices r ~ k" 1 play an important 
role in Eq. (94) and in the sum of Eq. (103). This results in the value of A kc (k), 
not in 2ir5(k), in the expression (95) for the power spectrum and (104) for the 
integral kernel, see also Eqs. (100)-(101). 

In order to improve the approach formulated above, let us replace, respectively, 
the Dirac delta-functions 2n5(k) and 2n5(k') with prelimit ones A kc (k) and 
A kc (k') in asymptotical Eqs. (123) and (124). With this, the expression for 
the kernel turns out to have the form, 



^{k,k) = 4 > cos (kr) — 

2 V, A c r 

r,r =1 L 

x sin(A;y) co y 
k c r' 

1-h 



[l-A fcc (A;)][l-A fcc (A;')]. (131) 

Let us compare the first expression (130) with exact Eq. (103) and zero- 
asymptotical Eq. (121). In the exact expression (106) for the kernel E(k,k') 
in the case of long-range correlations, fc c C 1, we admit that the following 
replacement 

I sw (k,k')^A kc (k)A kc (k') (132) 
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is correct. Note that such a model conserves all general properties of the func- 
tion I sw (k,k') described above after Eqs. (107)-(108), and as a consequence, 
provides all properties of the kernel E(k,k'), as well as the exact equations 
(31) and (113). 

The assumption (132) allows one straightforwardly to obtain the solution of 
the equation (113), 



2h 
1 - h 



F(k) = (2 + B sw ) 



1 + ^A^AOpA^A;) 



-(2 — B sw ) 



1 + 



1-h 



h 



A fcc (A:) 



h 



Now one should take into account the equalities 



1 + 



1 - h 



h 



A fec (A;) 



-i 



k c (l - h) 
k c h + 7r(l - h) 



A fec (A;), 



i + YA fc w]\w 



k ch+ %- h) ^^ 2 ^=i^k). 



Then the solution gets the form 



(133) 



(134) 



(135) 



2h 



T{k) = -(2 - B sw ) 



1-h 

+ 4k c h - (2 - B sw )[k c h - k c (l - h)\ 



k c h + tt(1 - h) 
From the normalization condition (42) one can obtain 

Ak r h 



2-5, 



(tt- k c )(l-h) + 2k c h' 



(136) 



(137) 



As a result, in the case of long-range correlations we have 



Hk) = 



2k c (l - h) 



(n - k c )(l - h) + 2k c h 



[A fec (A;)-l], 



F(r) 



2k c (l - h) 



(it - k c )(l - h) + 2k c h 



sin(/c c r) 
k c r 



(138) 
(139) 
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Evidently, this expression can be simplified by making use of a smallness of 
k c . Thus, the final correct expression gets the form, 



2k c 



[AjJAO-1], 



(140) 



7T 




2k c sm(k c r) 



-5„ 



(141) 



7r k c r 



V,o 



Note that the obtained result (140)-(141) is independent of h, hence it is also 
valid for h — 0. 



12 Concluding Remarks 

In this paper we suggest a method of construction of a binary Markov chain 
with a prescribed pair correlator. It is based on the integral relation between 
the pair correlator and the so-called memory function. The knowledge of the 
latter allows one to create binary sequences that have given short or long-range 
pair correlations. Our interest in this method is mainly due to its application to 
the problem of selective transport in one-dimensional disordered systems. As 
is known, the theory of quantum transport through finite samples described by 
a white-noise potential is fully developed. The key ingredient is the well known 
Anderson localization, that states that all eigenstates in infinite samples are 
exponentially localized regardless of the strength of disorder. The knowledge 
of the localization length of these eigenstates allows for a complete statistical 
description of all transport properties of finite samples, due to the single- 
parameter scaling conjecture. 

In contrast with the white noise disorder, the role of long-range correlations in 
random potentials is far from being properly understood. Recent results reveal 
quite unexpected properties of the transport in systems with such potentials. 
In particular, in Refs. [1,2,3,4,5,6,7,8,9] it was shown that specific long-range 
correlations give rise in a very sharp change of the transmission coefficient, 
when the energy of an incident wave crosses some value E c . This effect is 
similar to the mobility edge effect in three-dimensional solid state models, ac- 
cording to which on one side of E c all eigenstates are localized, in contrast 
with other side where the eigenstates are extended. It is assumed that analo- 
gous effects in one-dimensional geometry may find various applications, such 
as a construction of materials with a highly selective transport or explanation 
of a selective conductivity in DNA molecules. 

However, as was recently understood, an extension of this theory to binary 
sequences suffers from unexpected problems. We hope that our study of the in- 



29 



tegral equation (11) relating the memory function and pair correlator can shed 
light on the problem of binary correlated disorder. Since there is no general 
theory for solving integral equations, we have decided to start with particular 
cases allowing for analytical solutions. In this line, apart from the study of 
general properties of the main equation (11), we have solved relatively simple, 
however, physically interesting cases of short-range, exponentially decaying, 
and exponential-cosine correlations. In all these cases we did not met seri- 
ous mathematical problems. On the other hand, a more general problem of 
the step-function (95) that is of specific interest in view of the mobility edge 
transition, we were able to solve this problem for the cases when the mobil- 
ity edge is close to either one or other energy band. We would like to stress 
that these two cases are also interesting from the view point of different ap- 
plications, giving close analytical expressions for the pair correlators. We have 
found, however, that it is not clear how to obtain the general solution of the 
problem. It is still not clear whether the mathematical difficulties of solving 
of the integral equation (11) are principal, or they can be solved in a different 
approach. Note, that here we discuss the analytical results only, leaving aside 
a possibility to solve the problem numerically. We plan to do this in our next 
studies. 

To conclude with, we would like to note that the problem of construction 
of sequences with a power-law decay of the pair correlator, is quite specific. 
In particular, one should take into account strong restrictions imposed by a 
rigorous theorem stating that the Fourier transform of any pair correlator 
is a non-negative function (this corresponds to the fact that the Lyapunov 
exponent of the considered processes can not be negative). This results in 
a quite unexpected conclusions. For example, there are no sequences giving 
rise to a monotonic power decay for the pair correlator, such as K(r) = C/r 1 , 
where C and 7 are some constants. In other words, one should always take into 
account the oscillations of a pair correlator, such as sin(/c c r). These oscillations 
play an important role providing the non-negative values of the Lyapunov 
exponents (see, for example, Ref. [1]). In the latter case their influence on 
physical observables may be neglected. Therefore, it is not clear whether one 
can speak about "generic cases" of the power decaying correlations. 
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